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I. INTRODUCTION 

Generalized parton distributions (GPDs) [[1], § are considered to be a promising tool to 
study the partonic structure of hadrons (see Refs. [0-|6| for recent reviews). Hard exclu- 
sive reactions, which can be described in the theoretical framework of GPDs, provide us 
opportunity to access the information on GPDs experimentally. 

The importance of GPDs was widely realized in connection with the possibility to study 
the total angular momentum of partons in the nucleon. This allows to address the funda- 
mental question how the total nucleon spin is made up from the contributions due to quarks 
and gluons. The famous Ji sum rule |2j relates the appropriate Mellin moments of quark 
and gluon GPDs 1 H q,9 (x,^,t) and E q,9 (x,l;,t) to the fractions J 9 ' 9 of the total angular mo- 
mentum {i.e. the sum of parton spin and orbital angular momentum) carried by quarks of 
the flavor q and gluons respectively: 




dx x [H q (x, £,t = 0) + E q (x, «£, t = 0)] = 2J q ; 



/ dx [H 9 (x^, t = 0) + E 9 (x,£, t = 0)] = 2J 9 . (1) 
Jo 

Constraining J 9 from the experiment would add an important lacking "piece" to the 
tricky "nucleon spin puzzle". This require the knowledge of gluon GPDs H 9 , E 9 as func- 
tions of x for fixed values of £ and t. These GPDs can be most preferably studied in the 
hard exclusive electroproduction of J PC = 1 mesons (p ,u,<f)). For these meson electro- 
production channels H 9 and E 9 make contributions at leading order in l/Q (Q 2 refers to 
the initial photon virtuality) and in a s . Hard exclusive electroproduction of mesons and 
particulary the exclusive electroproduction of p° is now in focus of intensive experimental 
investigations (see e.g. ]7)). The important set of data have been already published [|J. 
Future experiments will provide even more precise data over a broader phase space. Thus, 
further theoretical development of GPD formalism that would allow the interpretation of 
this new experimental information is highly demanded. 

Unfortunately, the problem of GPDs extraction from the data is complicated by the 
fact that GPDs depend on several variables (longitudinal momentum fraction of partons x, 
skewness parameter £, momentum transfer squared t and the factorization scale). More- 

1 For the definition of nucleon GPDs we use the set of conventions employed in B. See below. 
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over, only the integral convolutions of GPDs with certain convolution kernels rather than 
GPD themselves enter the observable quantities. Thus, in order to extract the informa- 
tion on GPDs from the data, one of necessity has to rely upon different phenomenologicaly 
motivated parametrizations of GPDs and ingenious fitting procedures for the observable 
quantities. One of the most popular parametrizations of GPDs is the famous Radyushkin 
double distribution Ansatz (RDDA) [[| employed in numerous phenomenological applica- 
tions. In particular, the specific version of RDDA was adopted for gluon GPDs in the 



nucleon [|L0| [LJJ. 

There is the growing confidence that the present inability to describe some aspects of 
the available experimental data on the hard exclusive processes may be due to incomplete 
or inexact way of modelling GPDs (see e.g. discussions in [|12|-|T4|1). A possible alternative 
way to parameterize GPDs consists in employing of the so-called dual parametrization of 



GPDs |fL5|| . In this paper we apply the dual parametrization approach, which was originally 
developed for quark GPDs, to the case of gluon GPDs in the nucleon. We also perform the 
generalization of Abel transform tomography method |T(J for the case of gluons. Finally, 
we discuss some aspects of the fitting strategy for the hard exclusive processes observables 
based on the dual parametrization of GPDs. 



II. BASIC DEFINITIONS 



Following the conventions accepted in the unpolarized gluon GPDs in nucleon H 9 , E 9 
are defined as the Fourier transform of the matrix element of the nonlocal gluon operator 
between the nucleon states according to: 
1 f dX 



AXxP-n 



P-n J 2tt 
1 



(N(p')\F+»(-\n/2)F u + (\n/2)\N(p)) 



2P ■n U{P+ Y } [H9{X > e ' t)U ^ + 2^ E9{X > e ' ^^ A * 
The polarized gluon GPDs in nucleon H 9 , E 9 are defined as 
-i C dX 



P-n / 2vr \ w>\ 



^1—U(P + |) [H 9 (x,UKYj 5 4 



An/2)F+(An/2)|iV(p)) 
1 



E 9 {x,t,t) lb n^ 
2rriN 



U(P 



A, 



In (0), ([|) we employ the standard notations: n is the light-cone direction 



(2) 



(3) 

= o, 



P ■ n = P + ), P = + p'), A = p' — p, t = A 2 , the skewness variable £ refers to 



A + = — 2£P + ; F alS = |e a/37<5 F 7< 5 is the dual gluon field strength. Throughout this paper we 
adopt the light-cone gauge A + = 0, so that the gauge link does not appear in the operators 
in definitions (|J), @. 

The definition (0) differs from other definitions encountered in the literature (see |4|]). 
For < x < 1: 

H 9 (x,^t)\ here = H 9 (x,£,t)\ [3] + H 9 (-x,£,t)\ m = x (H 9 (x,£,t)\ [17] - H 9 (-x,Z,t)\ [17] ) . 
H 9 (x,^)\ h ^ e = H 9 (x,^)\ m - H 9 (-x,C,t)\ [3] = x (H 9 (x,£,t)\ [17] + H 9 (-x,£,t)\ [17] ) • 

The same relations hold for E 9 and E 9 respectively. 

As the gluon itself is its own antiparticle gluon GPDs H 9 (x,l;,t), E 9 (x,^,t) defined in 
(§) are even functions of x: 

H 9 (x, £, t) = H 9 (-x, £, t) ; E 9 (x, £, t) = E 9 (-x, £, t) . 

(4) 

Gluon GPDs H 9 (x,£,,t), E 9 (x,£,,t) defined in (|3|) are odd functions of x: 

H 9 (x, £, t) = -H 9 (-x, £, t) ; E 9 (x, t) = -E 9 (-x, £, t) . 

(5) 

Let us stress that in what follows we consider the gluon GPDs in nucleon (g), (|3|) on the 
interval < x < 1. 

In the froward limit gluon GPDs H 9 and H 9 reduce to usual forward gluon distributions 
in the nucleon g{x) and Ag{x), while GPD E 9 and E 9 are reduced to unknown gluon 
distributions, which we denote as e 9 (x) and Ae 9 (x): 

H 9 (x, 0, 0) = xg(x); E 9 (x, 0, 0) = xe 9 (x) ; 

H 9 {x,0,0) =xAg(x); E 9 (x, 0, 0) = xAe 9 (x) . (6) 

Note that the forward gluon distributions g(x), Ag(x) and e 9 (x), Ae 9 (x) are continued to 
the negative value of their argument according to: 

9(x) = -g(-x); e 9 (x) = -e 9 (-x) ; 

Ag(x) = Ag(-x) ; Ae 9 (x) = Ae 9 (-x) . (7) 
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The Mellin moments in the momentum fraction x are of major importance in the GPD 
approach. According to the polynomiality property of GPDs the x Mellin moments of gluon 
GPDs defined in (fj), ([3]) are polynomials of £. The coefficients of these polynomials are 
related to the form factors of the local twist two gluon operators: 

= S(-i)F tia i1^ tl K..i1^ t " n F^ . (8) 

Here, as usual, D* 1 is the covariant derivative, S denotes symmetrization in all uncontracted 
Lorentz indices and subtraction of the appropriate traces. More technically, the polynomially 
property for unpolarized gluon GPDs means that for odd N: 

„! N+ mod (AT,2) 

/ dxx> T - 1 H'(x,t,t)= V eh 9 Nk (t); 

J fc =o 

even 

-1 N+ mod (N,2) 

/ dxx N - 1 E^(x } ^t)= V ee 9 Nk (t) 
Jo k=0 

even 

with e 9 N>N+1 (t) = -h 9 NN+1 (t). (9) 
For the case of polarized gluon GPDs for even N 

/ dxx N - 1 H°(x,Z,t) = J2t k h 9 Nk (t); 

J fc =o 

even 

/ dxx N - 1 E*(x,U) = J2^Nk(t)- (10) 

J k=o 



ILL THE DUAL PARAMETRIZATION FOR GLUON GPDS 

Historically the first non-trivial phenomenological para-metrization of GPDs was the 
famous Radyushkin double distribution Ansatz (RDDA) suggested within the double dis- 
tribution representation of GPDs 0. In particular, a version of RDDA was adopted for the 
case of gluon GPD H 9 in the nucleon (@) 

H»(x, t, t = 0) = H 9 DD (x, + 29{i - \x\)iD<> {^j , 

(11) 
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where D 9 stands for the gluon D-term || and H g DD is built as a one dimensional section of 
a two-variable double distribution: 

H 9 DD (x,£)= / dp / da {5(x-/3-aO -5{x + p - a£)} h {b) {/3,a)pg{/3). 

JO J -1+0 

(12) 

The profile function h^ b '(/3,a) is parameterized as 

The parameter 6 characterizes the strength of £ dependence of the resulting GPD. The 
usual choice for the gluon case is b = 2. This is motivated by the interpretation of the 
a dependence like a meson distribution amplitude for hard exclusive processes. The cases 
6 = 2 correspond to the asymptotic behavior of a gluon distribution amplitude ~ (1 — a 2 ) 2 . 

The achieved theoretical understanding of both theoretical and experimental aspects of 
GPD physics hints (see e.g. |12, [13|, |18| ) at the necessity to introduce new GPD parametriza- 



tions which should be more general and flexible than the basic form of the RDDA employed 
in present-day mainstream phenomenology 

The possible alternative way to parameterize GPDs is the so-called dual parametriza- 
tion [15) (see 19— 22| for the recent development and discussion). Originally, the dual 



parametrization was formulated for the case of quark GPDs. Now we are going to generalize 
this approach for the gluon case. 

In the framework of the dual parametrization GPDs are presented as infinite sums of 
the t-channel 2 Regge exchanges Let us stress that the term "dual" is intended to lay 



emphasis on the natural association with the old idea of duality in hadron-hadron low energy 
scattering. The essence of the duality hypothesis for binary scattering amplitudes |3, |25 



can be summarized as the assumption that the infinite sum over only the cross-channel 
Regge exchanges may provide the complete description of the whole scattering amplitude in 
a certain domain of kinematical variables. 

More technically, in the dual parametrization the t-channel matrix element of the partic- 
ular non-local light ray operator O between the hadron states which enters the definition of 



2 The i-channel refers to the i-channel of the hard exclusive electroproduction reaction in question. E.< 
for the case of DVCS this is hadron pair production 7*7 — > NN . 
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GPD (see e.g. (Q), (|3|)) is presented in the following form: 

(N(p')N(-p)\6\0)~J2 E 7-175- x >W) iV (-P)l^> < WO|0> 



t — Mn 

f? 7 polarization _ _ _ . _ . . _ 

J of Rj RjNN effective vertex F.T. of DA of Rj 



(14) 



The sum in (|14|) stands over all possible t-channel meson resonance exchanges with suitable 
quantum numbers and of arbitrary high spin J and mass Mrj. For the classification of 
RjNN vertices see e.g. ref. |p6| . The distribution amplitude (DA) of the spin- J t-channel 
resonance occurring in (Rj \ O |0) matrix element is expanded in the eigenfunctions of the 
ERBL (Efremeov, Radyushkin, Brodsky, Lepage) evolution equation. For the gluon case 



these are the Gegenbauer polynomials Cfi (z) p7| . The on shell spin sum of spin- J t-channel 



resonance resulting from the sum over polarization in ( |14|) is expanded in the t-channel 
partial waves. 

Thus, in the framework of the dual parametrization GPDs have the form of double ex- 
pansions in the eigenfunctions of the ERBL kernel and in the t-channel partial waves. These 
formal series are then analytically continued to the physical region by means of the cunningly 
organized summation procedure. It is worthy to mention that the dual parametrization of 
GPDs shares many common features with the expansion of GPDs in collinear conformal 
partial waves |28| known in different versions [0, [14|, |29H3^1. The expansion of GPDs in 



collinear conformal partial waves arises naturally from the solution of the evolution equations 
to the leading order accuracy. In both cases the angular momentum of SO (3) partial waves 
expansion in the t-channel partial waves as well as the conformal spin appear as labels. 

The important point is that one has to take properly into account the complication 
introduced by the fact that the nucleon has spin-| . In order to be able to write down the 
formal series for the gluon GPDs in the nucleon in the framework of the dual parametrization 
it is necessary to point out the combinations of gluon GPDs H 9 , H 9 and E 9 , E 9 suitable 
for the partial wave expansion in the t-channel partial waves. This is an easy task since the 
Lorentz structure of the Fourier transform of nucleon matrix element of gluon light cone 
operators which enter the definitions (|2|), (^) is obviously the same as that of the quark 
light cone operator familiar from the definition of unpolarized and polarized quark GPDs 
respectively. Thus, in the complete analogy the case unpolarized quark GPDs [§, |33[ the 
partial wave expansion in the t-channel partial waves can be written for the electric and 
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magnetic combinations of unpolarized gluon GPDs 

H^ E \x, f , t) = H 9 (x, t) + T E°(x, C, t); (r = ' 



Am 2 N 



H 9 M(x, t ) = H 9 (x, £, t) + E 9 (x, £, t) . (15) 
and for the following combinations of polarized gluon GPDs: 

H 9 ^ (x, f , t) = H 9 (x, C, t) + rE 9 (x, £, t) . (16) 
One may come to the same conclusion employing the general method for determining 



the invariant amplitudes of a binary scattering process suggested in |35| (see also P5]] ). 
Using this method it is also straightforward to check that the electric combination H 9 ^ E ' 
and pseudoscalar combination H^iPS) are ^ Q ^ e expanded in the Legendre polynomials of 
cos Ot Pi{cosOt) while the magnetic combination H 9 ^ and H 9 are to be expanded in the 
derivatives of the Legendre polynomials P{ (cos 9t). Note that to the leading order of 1/Q 



expansion the i-channel scattering angle 3 is expressed through the kinematical variables as 



cosfl,= , +Q ). (17) 



1 

For the electric combination of gluon GPDs H 9 ^ J PC = ++ , 2 ++ , ... intermedi- 
ate meson states (e.g. /o, f%) contribute, while for the magnetic combination H a( - M ^ 
the contributions of J PC = 2 ++ , 4 ++ , ... intermediate meson states are relevant. For 
H 9 J PC = 1 ++ , 3 ++ , ... intermediate meson states contribute and, finally, for H 9 ^ PS ^ 
J PC = h , 2 h , ... meson states contribute. 

Taking into account all these considerations one can write down the following partial 
wave expansion for electric and magnetic combinations of unpolarized gluon GPDs in the 
nucleon 4 : 

oo n+1 / 2 \ / 2 \ ^ 



' / 2\ / 2\ z = / \ /i 



odd even 



(18) 



3 9 t is defined as the scattering angle in the center of mass frame of the i-channel of the hard exclusive 
electroproduction reaction (e.g. 77* — > NN for the DVCS). 

4 Note, that contrary to the case of quark GPDs in nucleon defined in (22| we do not introduce the factor 
"2" in the partial wave expansions for gluon GPDs. 
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00 n_ '~ 1 / 2\ / 2 \ 2 r: / \ /l ' 

H»<->(»,«,t) ^Etf'c^ 1 - f.) (i - 1) (f) * (i 



The partial wave expansion for the polarized gluon GPDs in the nucleon reads 

*'(^.')-Eg^)<i-S(i-S^-,(f)^(i 



*«w.o-gE«r'w(i-f)(i-?) , ^(f)flG 



(19) 



(20) 



n=2 Z=0 
even even 



(21) 

In fact, the most non-committal way to understand the partial wave expansions ([15), 



(PH), (p0|), (pl|) is to consider them just as formal series which satisfy the fundamental 
polynomiality property of gluon GPDs (|5|), ([TP]) . The summation procedure |T5| allowing to 
convert these formal series into rigorously defined expressions is reviewed in the Appendix [A|. 

In what follows we are going to consider in details the properties of unpolarized gluon 
GPDs in the framework of the dual parametrization. The summary of results for the case 
of polarized gluon GPDs in the nucleon is presented in the Appendix |B[ 

Let us first consider the electric combination. For odd iV the (N — l)th Mellin moment 
is indeed the polynomial of £ of order N + l 

1 N+l 

dxx»- 1 H'W(x,t,t) = £?h'W(t) 



k=0 
even 



P N V V R* <*> (tV P (1\ U{l+n) {2 + H) {3 + H) T M^1 (22) 

it it \U 9-2^r(i + ^)r(| + ^f±^) " { ' 



n =i 1=0 v- 1 - 1 2 / v2 

odd even 



The set of coefficients h 9 J E \t) can be expressed through the generalized form factors B^ E ^ 
according to 

N n+l r {2-k+l+N^ 
k+l-N-1 



u N,k \ L ) Z^Z^^nl \ L )\ L ) g nk+AY( l+k + l ~ N ) T(2 fe | A r ) 

odd even 

7i (1 + n) (2+n) (3 + n) r(JV) 

X r (2 = ^fiV )r( 5+^ • 



(23) 
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For the magnetic combination the (N — l)th Mellin moment (JV- odd) is the the polyno- 
mial of £ of order iV — 1 

-1 iV-l 

/ ^-^W(x,e,t) = ^« ( f } W = 

even 

^ V V h« <") r/v p' f M w(1 + w) (2 + n) (3 + n) EffiEW r ^ 
t^tr u/ 9-2^r(i + ^)r(| + ^±£±iV) • 



n=l l=o \ ^> s " " *■ V * 1 2 / V 2 

odd even 



The corresponding set of coefficients h^ k (t) is expressed through the generalized form 



factors as follows 



N n+1 t_ 1 + k _ N ) T ( 2-k+l+N ^ 



"iv,fc w 1^2-^ nX y,y ' o 2 < ^ +1 r( l+fc +^ Ar )r(2 fc i A r ) 

odd even 

n (1 + n) (2 + n) (3 + n) r(JV) 



r(^v)r(5±^) 

(25) 

Gluon GPDs H 9 ( E,M >(x,£,t) are normalized according to: 

fdxH^ E \x^t) = Mi(t) + ~(1 -rK(*K 2 ; 

/ dxH 9(M) (x,£,t) =2J 9 {t). (26) 

Here Mf (i) stands for the t-dependent momentum fraction carried by gluons in the nucleon; 
J 9 (t) denotes the t-dependent fraction of angular momentum carried by gluons and d 9 (t) is 
the first coefficient of Gegenbauer expansion of the gluon D-term (p^). 

IV. THE PROPERTIES OF H^ E ' M \x,£,t) IN THE DUAL PARAMETRIZATION 

To sum up the formal series for the electric and magnetic combinations of gluon GPDs 
(fL8|), (0) we employ the techniques developed in |l5j (see also discussion in f22|). Some of 
the additional technical details specific for the gluon case are presented in the Appendix [A|. 

To proceed with the summation of the formal series (|18D, (|19|) we introduce two sets of 
gluon forward-like functions G 9 2 l and Gfi , whose Mellin moments generate electric and 
magnetic generalized form-factors B^ E ' M \t): 

B^Zlit) = f dxx"G 9 ^ M \x, t) . (27) 
Jo 
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The explicit expression for the electric combination of gluon GPDs H g ( E > through the 
corresponding forward like functions reads: 

00 p2v 
v=0 



(28) 



Note, that the second term in (|28|) is the pure D-term contribution. The result for the 
magnetic combination of gluon GPDs H 9 ^ can be obtained in the similar way applying 
the differential operator ^1 — x-^ — £jQ : 



00 p 

E 

u=0 



i 1 - x i " t" [H9 {M) M (a; ' ^ t] + H9 m M e ' t)] 



(29) 



The functions H9( E > M )W(x,£,t) which appear in (p8|) and (|29|) defined for — £ < a; < 1 are 
given by the following integral transformations: 



6{x > 0- / dy 
71 Jyo 



9 1 2d 2 



+6(\x\<0- / dy 

Jo 

2 \ 2 2j/-3 



9 1 2 9 2 



.r 



e»v e) ^ Cl - 1 - 3 U J m [u (2u-i)(2u-i+i) 



G% M \y,t) 
G% M \y,t) 

6y 2v-l-2 



*2 



*1 





S 2 ) 


^2x s - x\ 


-e 


J x 2 s - 2v (l 




d.M 2 L_ 





, ^2 Xs -xi-e 



(30) 



with P_ n (x) = P n -i(x)- Here, as usual, x s = 2jj 



x—£s 



Si, U 



1, ...4) stand for the four 



roots of the equation 2x s — x\ — £ 2 = (see ( |A12| ) for the definitions) and y is defined in 
( |A14| ). The integrals in ( |30"D are well convergent for the set of the forward like functions 



with the following small- y behavior G2 U (y) ~ ysu+s, with a < 2. 
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Point £ = x 



The important limiting case in which the expressions (|28|), ( p9[) are reduced to much 
simpler forms is the point x — £. 

u=0 



H° <*> & 1 1) = ±Z f *L jr y*Gg> (y, t) 

J 1 y u=0 



f - ^ - 1 



1 



(31) 



Point £ = 1 

Another notable limiting case is £ = 1: 



H 9{E \x,£= l,t) = -(1 - ./■ 



2 ) 2 E / d w 



(l — 2xy + y 2 )z (1 + 2xy + y 



2^ 



2i/-2 



;r 



even 



G£ 3 (j/,t); 



(32) 



and 



^W(x,£ = l,t) = ^(l-x 2 ) 2 E / ^ 

,,— n •'O 



i/=0 



1 



(1 - 2xy + y 2 )^ (1 + 2a;?/ + y 2 )? 



J=0 
even 



y dy [ y dy G ^ {V ' t} 



(33) 



Forward limit and G^ E,M \x,t 



Let us introduce the convenient notations for the combinations of the t-dependent parton 
densities to which GPDs H 9 ^ E ' M ^ are reduced in the limit £ — > 0: 

H 9{E \x,£ = 0,t) = xg{x,t) + Txe g {x,t) = xg (E) {x,t) ; 
H 9(M \x^ = 0,t) = xg(x,t) + xe 9 (x,t) =xg^ M \x,t), 

(34) 
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where g(x,t) and e 9 (x, t) stands for the t-dependent gluon distributions in the nucleon 
(xe 9 (x,t) = E 9 {x,^ = 0,t)). Employing the general results (p3|), (|25| ) one can check that for 
odd N: 

j\xx^H 9 ^{x^ = 0,t) = j\xx N g^(x,t) = {N ^n + 3) ^ ] 

J\ xx ^H^i = o,t) = [dr,\r^,,n = {N + ^ N + *l N + 3) <^(t) 

(35) 

Inverting the Mellin moments in ( pop we express the forward like functions Gq' related 
to the t-dependent gluon distributions 

d V JE)/„. +\ o_ f n (E)t t \ 



Gf\x,t) = 9x* [ *y g W(y,t)-3x ! %g ( 

j x y j x y 



: (M), ,v = _9 2 f 1 dy, M)(ni ,x ^ ~ f 1 ^JM), ,x , 3 f 1 dy m 



G^\x,t) = --x' ^ M) (y,t) + 3x ^(y,f) + J ^g (M) (y,t). 
2 Jx y s Jx y 2 J x y 



(36) 



The normalization is 



J dxxG$ S) (x,t) = ~M°(t); 

/ dxxG ( M \x,t) = ^2J 9 {t). (37) 
Jo 8 

Thus the information on the fraction of nucleon total angular momentum carried by gluons 
is encoded in the magnetic forward like function Gq . 



Small £ expansion 

It is helpful to consider the expansion of electric and magnetic combinations of nucleon 
gluon GPDs in powers of £ around £ = for fixed x (x > £). For the electric combination 
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jjg(E) the corresponding expansion to the order £ 2 it is given by: 



1 

32 



~Tx + + ^ (1 + ^^(M) + ^ - ^T^HM) 



128 \y 



x 



+ 



1 

24' 
3 



1 

24" 



d 2 



<9x 2 



128 



3 1 
+ 



y 2 \ 128 \x 



128 \y 



x 



dyGf\y,t) 



35 (yyk 5 



x\ 5 3 
128 \y) + 256 Vj/ 



3\ "1 

X N 2 



0(£ 4 ) ■ 



256 Vx, 

For the magnetic combination H 9 ^ the expansion to the order £ 2 reads 

H' {M) (x,i,t) = \xGi M \x,t) - l x ^G' ' J \x,t) + ^C'"°(i,') 



^4"W) + ^(i-5**)|;Grw) 



5 

48 



15 fx 
256 

_L 2 — 

48 X dx 



15 / x 



y \2 15 /'ac 



35 /y\ 2 15 /x 



15 /x\ 2 



512 Vy 



256 \y 



+ o(e 4 )- 



(38) 



(39) 



Modelling H 9( - E \x,£,t = 0) 



On fig [T] we show the results of the numerical computation of G contribution into the 
electric combination of gluon GPDs H 9<yE ^ at t = for several different values of £. As 
the numerical input for forward gluon distributions we used MRST LO fit for nucleon 
parton distributions at Q 2 = 1 GeV 2 . 
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0.2 0.4 0.6 0.8 1 

X 



FIG. 1: Electric combinations of nucleon gluon GPDs H 9 ^ E \ at t = for different values of £. As 
the numerical input for forward gluon distributions we use MRST LO fit |}(| for nucleon parton 
distributions at Q 2 = 1 GeV 2 . 



V. ELEMENTARY GLUON AMPLITUDES AND ABEL TRANSFORM TOMOG- 
RAPHY 



The typical convolution integral involving gluon GPDs F 9 = {H 9 , E 9 } relevant for the 
calculation of hard exclusive J PC = 1 meson electroproduction at the leading order in 
1/Q and a s reads [|J : 

1 1 



Jo x 



£ — x — ie ^ + x — ie 



(40) 



We define electric and magnetic elementary amplitudes 



x 



£ — x — ie £ + £ — ie 



(41) 



In order to obtain the partial wave expansions of the elementary amplitudes in the t-channel 
partial wave one has to substitute the formal series (0), fll9D into the corresponding convo- 
lution integrals. For odd n > 1 



x 



2\ 2 



<M 1 " £2 J C n-1 1 7 ) Z 



1\ L 



£ — x — ie £ + x — ie 



4 
3 



(42) 
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Thus 



A 00 n+1 /I \ 

n=l 1=0 \ > / 

odd even 

^ (M, «.t) = iEE B 'l M, M?) ■ 

n=l i=0 ^ \> / 



(43) 



These formal series can be summed exactly as for the case of singlet electric and magnetic 
quark GPDs. The form of the resulting expression actually differs only by a factor |. The 
expressions for gluon electric and magnetic elementary amplitudes read 



6J o X „=o 



l-^ + x 2 



+ 



-25, 



2v n g(M) 
2v 1 



x,t) 



1 + f + X 2 



+ 



(44) 



1 - f + x 2 



-25, 



1 + f + x 2 



(45) 



We introduce electric and magnetic gluon GPD quint-essence functions: 

00 

N^ E \x,t) = J2^ G{°\x,t); 



u=0 

00 



(46) 



u=0 



The imaginary parts of gluon electric and magnetic elementary amplitudes then read: 



Iim4 9(M) (£,t) 



dx 



3 / i-Vi-e 2 x 
£ 



N^ E \x,t) 



f-x 2 -l 



(47) 



3 liWHE 

£ 



dx 



d N^ M \x,t) 
dx 1 — £x 



f - x 2 - 1 



(48) 
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The real part of the elementary electric gluon amplitude is given by 



3 Jo x 



1 _2x +x 2 



1 + f + x 2 



x z 



dx 



+- / — N g W(x,t) 



l + ^ + X> 



2(1 - r)D 9 {t) , (49) 



where D 9 (t) stands for the gluon D- form factor (|5l|). Finally, the real part of the magnetic 
gluon amplitude reads: 



2 r — i — 



dxy/l + x 2 



1 



+ 



1 



1 + x 2 - f Jl + x 2 + f VI + x 2 



X ■ 



d VT+x 1 



dx \ 1 — x 2 



N 9 ^ M) (x,t) 



3 J i-VT^ 

s 



dxVl + x 2 



l + x 2_2x VT+ 



x A 



d (\fTTx 2 
dx \ 1 — x 2 



N 9 ^ M \x,t) 



The gluon D-form factor is defined according to 



D 9 (t) = Y,d 9 n (t) 



n=l 
odd 



-1 



dz- ^D 9 (z,t) 

1 — z A 



(50) 



(51) 



where D 9 (z,t) stand for the gluon D-term for which the following Gegenbauer expansion 
was adopted in ||: 



D 9 (z,t) = -(1- : 



2\2 



d\(t) + di(t)Ci(z) + dl(t)Cl(z) + 



(52) 



According to the analysis presented in [[13], |37|, the real part of the elementary amplitude 
A 9 ^ can be expressed through its imaginary part with the help a dispersion relation with 



one subtraction in the variable oj = | for the fixed value of t: 



A 9 ^(U) = 2D 9 (t) + - f\' (— 1 

k Jo \?-? 



ImA 9 ( E \£ -ie,t) 



(53) 
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The subtraction constant is given by the value of the amplitude at the non-physical point 
u = (C, = oo). It is known to be fixed by the D-term and equals 2D 9 (t). 



In |T2], [13], [32| it was suggested to fix the subtraction constant in the dispersion relation 
(^3]) in terms of the imaginary part of the corresponding elementary amplitude assuming 
special analytical properties in j of the function 

d>(j) = 2u (t) = / dxx*- [H^ E \x,x,t) -H^ E \x,Q,t)] . (54) 

For entire odd j > 1 ( |54"D defines in the family of sum rules for the specific combinations 
of coefficients ^i+j 2u a ^ powers of £ of the Mellin moments of GPD H 9 ^ (see eq. ( p2|) 
for the definition). For j = —1 the integral is divergent since it is usually assumed that 
j-(H 9( - E \x,x,t) — H g ( E \x,0,t)) ~ l/x a with a ~ 1. In order to provide the desired ex- 
pression for the D- form factor fl54|) is to be properly analytically continued to j = — 1. To 
make it explicitly, let us assume that -(H 9 ^ E '(x, x, t) — H 9 ^ E '(x, 0, if:)) belong to the class of 
functions with power like behavior for x ~ 0, which can be presented as the finite sums of 
singular terms PU| : 

F: F(x) = J2-^fr(*)- (55) 



r=l 



We suppose that for all r = 1, ...R a r < 2 and f r (x) are arbitrary functions of x infinitely 
differentiable in the vicinity of x = 0. It is also supposed that f r {x) have zeroes of a 
sufficiently high order for x — 1. Then the relation for the gluon .D form factor reads as in 

oo „! 

(56) 

The lower integration limit "(0)" in ([56]) symbolize that we use the so-called analytic (or 
canonical) regularization | |40| |: 

f dx 4^ = f dx -L [f( x) _ /(o) - x f(0)] - /( )i - /'(0)^- (for a < 2) . 

(57) 

In the framework of the dual parametrization assuming that N 9 ^(x,t) and G (E \x,t) 



also belong to class (^) (see discussion in [Q) one may rewrite the sum rule fl56"|) for the 
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gluon Z)-form factor as 



D»(t) 



2 r 1 dx 



-Gf ) (x,t)(-^=-l)+^ / - N°M(x,t)-G^(x,t) 





r 1 dx - 








1(0) % . 



3J ^ " VVlT^ J 3 J, 0) x L ^ ' ' UV "J yiTa; 



2 

(58) 



Thus the maximal amount of information on gluon GPDs which can be extracted from 
the leading order gluon amplitude of hard exclusive electroproduction of J PC = 1 mesons 
can be quantified in terms of gluon GPD quintessence functions N 9 ^ E,M ^ and the value of 
gluon D- form factor D 9 or, equivalently, in terms of N 9<yE,M ^ and Gq . The method of the 
Abel transform tomography for the dual parametrization of GPDs suggested in |Tj| allows 



to invert the integral convolutions for the imaginary parts of the elementary amplitudes in 
order to recover GPD quintessence functions. The generalization of this method for the 
case of gluon GPDs is straightforward. The expressions for gluon electric and magnetic 
GPD quintessence functions through the imaginary parts of the corresponding elementary 
amplitudes (|4lD read: 

N 9{E \x,t) = 



2 H 1 + X 2),J_^ {i -— 2 

3 + (E) (J2x 

2vr v^T^ \1 + x 



2 - 




It is also extremely instructive to consider the small-£ asymptotic behavior of 
lmA 9 ( E ' M \£) in the framework of the dual parametrization. Assuming the power law behav- 
ior 

i x ) ~ — for gluon electric and magnetic GPD quintessence functions for small x 
we obtain for £ ~ 0: 



ImA 9(E) (£) 



2"+i r(i)r(a + i; 

e 3r(« + i) 



W^(£) „ g^£(i)£(^+i) . (61) 
£« 3T(a) V ; 

The case of particular importance is the so-called minimalist dual model in which only 
the contributions of the forward like functions Gq' M \x) ( j36[) are taken into account. Let 
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us assume the power law small- x asymptotic behavior for the electric and magnetic gluon 
densities 

^(*)~JL. (62) 
For the forward like functions Gq (x) the power law behavior ( |52] ) results in 



3(1 + 2a) 1 
(a + l)(a + 2)aj° ' 
„.!/,, . 3(1 + 2a) 1 



Gq (x) ~ . y —, ' oX — . (63) 
u v y a(a + l)(a + 2) x a y ' 

The contributions of the forward like functions Gq (x) into ImA 9 ( E > M ) have the following 
asymptotic behavior for £ ~ 0: 

I m ^ M 'K)- 2 ° + ^ )r( \ +a) f. (M) 
Go vs; r(3 + a) £ Q v ; 

VI. SKEWNESS EFFECT FOR SMALL £ IN THE DUAL PARAMETRIZATION 

In this section we discuss some aspects of GPD modelling in the framework of the dual 
parametrization. The important characteristics of the particular GPD model is the so-called 
skewness effect (also known as the skewness ratio). In JL3|, [14| it was suggested to define the 
characteristics of the skewness effect for quark and gluon GPDs H q ' 9 for t = as ratios of 
the GPDs on the cross-over line x = £ to the appropriate parton distribution functions, to 
which GPDs are reduced in the forward limit. Thus, for quark GPD H q the skewness ratio 
is defined according to 

q _ H«(x = Z,Z,t = 0) _ #"(£,£, 0) 
' " tfff(£,0,* = 0) g(£) ' 

(65) 

while for gluon GPD H 9 it reads 

g H°(x = Z,Z,t = 0) _ H°(Z,Z,0) 
T H*(£,0,t = 0) £g(0 

(66) 

It is usually assumed that the small-£ the asymptotic behavior of quark GPD on the 
cross over trajectory and the corresponding PDF is governed by the same leading Regge 
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trajectory: 



H*(Z,£,t = 0) = -lmA q ^,t = 0), q(£) 



(67) 



Similarly, for the gluon case: 



^,e,t = 0) = ilm^t = 0) ; «?(£) 

(68) 

For £ ~ this makes skewness ratios (|55|), fl6"6"D independent of £. The physical value of 
a q ' 9 (t = 0) in the relevant kinematical region is given by: a q,9 [t = 0, Q 2 = 4 GeV 2 ) = 
a 9,9 = 1.1 -T- 1.2. Obviously, the model without skewness in which H q (x,^) = q(x) and 
H 9 (x,l;) = xg(x) corresponds to r q ' 9 = 1. 

Quark skewness ratio can be directly related to the observable quantities. Since this point 
for some time was a rather knotty problem in the literature we present certain details of the 
relevant calculations following Ref. fl4] . It is worth to mention that often in the literature 



somewhat different quantities are employed to indicate skewness effect at small £. They are 
defined as the following ratios [|J: 

~ 9(20 ' H ~ 2^(20 • tb9j 

The relation between the two definitions of skewness effect under the assumption of Regge 
like asymptotic behavior fl67"D , (|68|) is given by 



R q = 2 a" r Q . R g = 2« 9 "V . ( 70 ) 

The definition (|6^) is inspired by the fact that the observable ratio of DVCS and DIS cross 
sections is reduced to R q rather than r q . Let us present some details concerning this issue. 



In Ref. HT[ in order to characterize the magnitude of skewness effects present in the DVCS 



process for small xg» the following ratio R was introduced as observable quantity: 



r(q 2 , w) — (71) 

where W stands for the total invariant energy, b(Q 2 ) is the fitted t-slope parameter 
(da Dvcs /dt ~ e~ 6 '*'). In the LO analysis the DVCS cross section a^vcs in the small 
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XBj regime is governed by quark exchange mechanism and is known to be dominated by the 
contribution due to the imaginary part of the amplitude: 



a2 e.m. nX Bj ,x A + _a\|2 1 , 2\ 



vdvcs(xbi ~ 0, Qr) ~ em Q4 a3 \ImA DVGS (Z, t = 0)| z ^y(l + p 



(72) 



The imaginary part of the leading order DVCS amplitude is given by Im A nvcsi^t = 0) = 
tiHdvcs(£,, 6, t = 0), where Hdvcs is the combination of the singlet unpolarized quark GPDs 
Hdvcs = ^H u +\H d . The parameter p 2 in (|n]) , ( [72]) refers to the small correction due to the 
real part of the amplitude. Since the virtual photon is assumed to be transversely polarized, 
in the case of DVCS it is also to be taken transversely polarized in the DIS amplitude in 
(0). The text book expression for the transversely polarized DIS cross section reads: 



<M7 V -+ x ) - ^ oi e . m . = 47ra e . m . — J - . (73) 



Now, since for small x B j the skewness parameter is given by £ ~ xbj/2, the definition (|7T|) 
for the observable R can be rewritten as 

Hdvcs (^,^,t = 0) 
H D vcs(xBj,t = 0,t = 0) 
Contrary to the original statement of Ref . [EI] , the model without skewness corresponds 
to R = 2 aq r q ~ 2. Thus, as it is emphasized in [|13| f[4| , the results of the experimental 
measurements of R performed by HI Collaboration [4l|| presented on Figure ^unambiguously 
allude to no skewness effect for small Xsj- 

Let us now discuss skewness effect in the framework of the dual parametrization. Em- 
ploying the results of the previous section and of || we compute the skewness effect for 
small £ in the minimalist dual model in which only the contributions of the forward like 
functions Q^f' M \x) and Gq M \x) assuming Regge like behavior (|67|), ( |68D of input electric 
and magnetic combinations of parton densities. The corresponding result reads: 



q(E,M) 



2 aq T(a q + -) 
~ prr \ ^- w 3/2 for a q « 1 ; (75) 

^ r(§)r(2 + o«) 1 > k * 



r 



g(E,M) 
G 



2 a9+1 T{a9 + |) 



1 for a 9 « 1 . (76) 



^ r(|)r(3 + a») 

It is extremely instructive to compare the skewness effect ( [T6|) in the minimalist dual 
model to that in the commonly used version of RDDA with the same asymptotic behavior 
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FIG. 2: The observable ratio R (fn]), shown as a function of Q 2 for fixed W = 82 GeV. This 



Figure is taken from [41] 



), (|68|) of the input parton distributions. The result for the corresponding skewness ratio 



presented e.g. in [T4| reads: 



1 RDDA 



RDDA 



H 9 DD 



2 2b - aq T(b+l)T(l + b-a q ) 
r(|)r(2 + 26-a«) ' 

2 2b+1 - a3 r(b + f )r(2 + b - a*) 

r(|)r(3 + 26-a9) 



(77) 



Comparing this to (7^) and (|76|) we conclude that both in the quark and gluon cases RDDA 
with b = a q ' 9 results in the very same skewness effect for small £ as the minimalist dual 
model with the same input. 

Thus, in the quark sector the minimalist dual model produces the non zero skewness 
effect for small £. It makes this model unpractical for the description of DVCS at small XBj 
due to systematical ~ 50% overshooting of the available HI cross section data. At this point 
the minimalist dual model has no advantage comparing to the common version of RDDA 



plagued by the same problem (see [j|2] and the discussion in [131, |14|) At the same time, in 
the gluon sector the skewness effect in the minimalist dual model turns to be very small. It 
is worth to mention is the growing confidence that the gluon skewness ratio is also not given 
by the conformal ratio ( ff6|) (see discussion in |I4]). The pure LO analysis of hard exclusive 
p° electroproduction rather hints that r 9 < 1. 

In fact one may recognize in (^) the well known result obtained in []38 using the Shuvaev 



integral transform approach G9|. In [14| it was suggested to call it as the conformal ratio 
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since this is nothing but a Clebsh-Gordan coefficient occurring in the conformal partial wave 
expansion. It is still unclear if the skewness effect ([73) is an indispensable feature of GPD 
phenomenology (see discussion in [JH, |38|). 



Let us emphasize that the skewness ratio in the framework of the dual parametrization 
is not with necessity given by (JT6|). In order to alter the skewness ratio one has to take into 
account the contribution of subsequent forward like functions Q2V {G2u) singular enough 
to change the small £ asymptotic behavior of if"(£,£). For this one has to assume that 
Q2u(%), G2v(x) ~ l/x 2u+aq ' 9 . The rigorous way to handle the occurring divergencies of the 
generalized form factors B 2u -io 



Bt E ^\t)= f dxx^Q^ M \x,t) 
Jo 



Bl^\t)= / dxx^G% M \x,t) (78) 
Jo 



and of the quark and gluon .D-form factors was described in details in . 

In the remaining part of this section we consider a toy model in order to briefly sketch 
the possible fitting strategy for hard exclusive scattering observables based on the dual 
parametrization of GPDs. For definiteness we discuss the case of DVCS at the leading 
order. 

In order to fit the experimental data we have to propose certain Ansatz for the relevant 
GPD quintessence functions N(x,t). The contribution of the forward like function Qo(x,t) 
into GPD quintessence function is entirely fixed in terms of the t-dependent parton dis- 
tributions. Since for the moment we discuss DVCS let q(x,t) be the DVCS combination 
of t-dependent singlet parton distributions q(x,t) = ^u + (x,t) + |<i + (x,t) . The observable 
quantities can be expressed in terms of the standard elementary amplitude 



A(£,t)= [ dxH(x,£,t) 
Jo 



£ — x — it £ + x — it 



(79) 



where H(x,£,t) is the DVCS combination of singlet quark GPDs. 

The leading singular behavior for x ~ of the t-dependent PDF q(x, t) is assumed to be 
determined by the linear Regge trajectory a{t) = a + a't: 

q(x,t)~ c ^, (80) 
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where c q > is the numerical constant. The value of the intercept in the relevant kinematical 
domain is a (Q 2 = 4GeV 2 ) = 1.1 1.2, The slope parameter a' can be fixed e.g. with the 
help of the form factor sum rule |39[]: a' = 1.1 GeV~ 2 . 



The alluring possibility is to take advantage of the opportunities provided by the Abel 
transform tomography method and instead of modelling GPD quintessence function employ 
a model for the imaginary part of the elementary amplitude lmA(C,, t). The real part of the 
elementary amplitude Rev4(£,t) can be then computed with the help of GPD quintessence 
restored from Imyl(^,t) using the tomography method. The important advantage of this 
approach is that the amplitude computed in this way possess proper analytic properties 
in £ and automatically satisfies the fixed t dispersion relation with one subtraction in the 
variable u = ~ |33], |37| : 

m t) = W( t ) + i [ < - ) IraAtf - it, t) , 

(81) 

since the expression for ReA(£, t) in the dual parametrization through the GPD quintessence 



function is equivalent to the the dispersion relation (81) . 

Let us introduce the following notations: lmA®° and lmA N ~®° for the contributions to 
the imaginary part originating from Qo(x, t) and N(x, t) — Q (x, t) respectively. Analogously 
for the D form factor we introduce the notations D®° and D N ~®° . 

The leading singular behavior of ImA^ ^, t) computed from the Qo{x,t) corresponding 
to t-dependent PDF q(x,t) with the asymptotic behavior (^) is 

Im^f£ t) ~ 2QW+1 + I) = cQo(t) J_ f82) 

[mA {U} c r( (t) + 2) - c (t) ew (82) 

where we define the function 

2«(*)+ir(i)r(a(*) + 1) 



C Qo {t) =c q 



T(a(t) + 2) 



C Qo (0) ~ yc 9 for awl. (83) 
The contribution of Qo to the D form factor can be computed using 

D*W = (^ - l) • (84) 

The skewness effect in the model which includes only the contribution of Qo is given by 
the conformal ratio ([75]) . In order to make skewness effect satisfy HI measurements 



we 
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need to tune the small £ asymptotic behavior of ImA(£,t = 0). For this issue let us require 
the following asymptotic behavior of lmA N ~®° for small £: 

1mA N -*°(£,t = 0) ~ - l -C^{t = 0)^ = £7^(t = 0)-L . (85) 

This choice makes the skewness effect consistent with HI results. 

In order to parameterize the effect of the contribution oi N(x, t)—Qo(x, t) to the imaginary 
part of the DVCS amplitude one may try to employ the following class of functions: 

l mA N-Qo^t) = tf"-«o(t)JL(l _£)/». (86) 

The corresponding truly non-forward contribution to the GPD quintessence function 
N(x,t) — Q Q (x,t) can be recovered form fl86|) employing the standard Abel transform to- 
mography procedure and used to compute the corresponding contribution to the real part 
of the elementary amplitude KeA N ~^°(^,t). 

Finally, the value of the N — Qo contribution to the D form factor can be computed 
rule employing the analyticity assumptions with the help of inverse momentum sum rule 

D N - Qo (t) = [ ^■lmA N - Qo (£,t) = C N - Qo {t)B{l + P,-a{t)) , (87) 

J(Q) £ 

where B is the Euler beta function. 

The analytical properties of D N ~®°(t) in the variable t require much attention. The 
Euler beta function in (|87[) has poles in t for t = ~ a J~ k with k = 0, 1, .... There is a finite 
number of "tachion" poles at negative values of t which do not match with the t-channel 
resonance exchange picture forming the basis for the dual parametrization approach. E.g. 
for 1 < a < 2 there are two "tachion" poles at t = — and t = — 9Lz r. This can be seen 

" a' a' 

as an indication that the simple Regge motivated substitution is inadequate for modelling 
the t-dependence of the DVCS amplitude. A more sophisticated Ansatz with the non-trivial 
interplay between the £ and t dependencies is needed for this issue. 

The interesting possibility is to employ the general form fl5E| ) assuming the special form 
of the t-dependence of C N ~®°(t) in order to get rid of "tachion" contributions into the D 
form factor: 

C N-Q {t) = C N-Q {Q) f t+ ^( t+ 2L=l) . (88) 

V a! J \ a' J a (a — 1) 
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This results in the following expression for the D form factor: 

v ; v J a{a- 1) T(l + /3-a(t)) v ; 

The value of C N ~®°(0) is necessarily negative in order to reduce the skewness effect for small 
£ (see (j85|)). It is interesting to note that for (3 > a(0) — 1 the value of the D N ~®° form 
factor computed from fl89|) turns to be negative for t = 5 . It is explained in ^3|, f|3| that 
the negative value of the D form factor is intimately related to the spontaneous breaking 
of chiral symmetry in QCD. Some arguments in favor of the negative sign of the D form 
factor are also presented in j44j in the framework of the simple model of very large nucleus. 



In this case the first coefficient of the Gegenbauer expansion of the D form factor can be 
related to the surface tension of the hadron medium and its negative sign follows from the 
requirement of the mechanical stability of the system. Thus, at this point our model turns to 
be consistent with the achieved theoretical understanding of the underlying physical picture. 
Finally it is interesting to note that asymptotic behavior of the D form factor D N ~®° (t) for 
— t — > oo is power like: 

D N - Qo (t) ~ (-t)-P +1 . (90) 

This simple exercise shows the power of the analyticity assumptions allowing to unambi- 
guity fix the subtraction constant in the dispersion relation ( |81"D in terms of the absorptive 
part of the amplitude. The value of the D form factor turns to be determined by the small £ 
asymptotic behavior of the DVCS cross section. It is extremely important to check whether 
this scenario is consistent with the available experimental data. 

VII. CONCLUSIONS 

In this paper we consider the application of the dual para-metrization approach to the 
case of gluon GPDs in a nucleon. We construct the partial wave expansion for both unpolar- 
ized and polarized gluon GPDs in the nucleon and present the explicit form of the integral 
transform allowing to rigorously sum up these formal series. We present the expressions 
for the elementary leading order amplitude entering the description of hard exclusive meson 



One may check that the contribution of Qo (84) into the D- form factor is also a small negative number. 
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production in the GPD formalism. We also discuss the generalization of Abel transform to- 
mography approach for the case of gluons. We argue that the dual parametrization provides 
opportunities for the more flexible modeling GPDs in a nucleon. We strongly suggest to 
use the fitting strategies based on the dual parametrization to extract the information on 
GPDs from the experimental data. Let us also stress that in principle there is no essential 
difference between the dual parametrization of GPDs and the Mellin-Barnes representation 



with the expansion of Gegenbauer moments in the in the t-channel SO (3) |T2|, [L4[ (see dis- 
cussion in 



45|). However, the full inversion formula relating the two parametrizations is still 



unknown. 
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A. SUMMING UP THE FORMAL SERIES FOR UNPOLARIZED GLUON GPDS 



The method suggested in [15] for the summation of the formal partial wave expansions of 
the type for GPDs in the framework of the dual parametrization consists in presenting 
GPD as the result of convolution of a certain convolution kernel with the set of forward like 
functions whose Mellin moment generate the generalized form factors B n [. 

The explicit form of the corresponding convolution kernel was presented in [13] for the 



case of quark GPD (see also |22j for the refined version of the derivation). In this appendix we 
present the summary of relations employed for the derivation of the integral transformation 
([28]) expressing GPDs H 9 ^ through the set of the forward like functions G^' . 
We introduce the common useful variable 

* = 2 <r^" (A1) 



28 



with < y < 1 and consider the discontinuity 6 

*1 1 „2 / i\JV-l /•! 



Employing the Rodriguez formula for the Gegenbauer polynomials Cn (z) ||46|| : 

(_l)n-i / 5n n 2\ / 9 



one can derive the following basic relation 

(, a 1 2 B 2 \ f 1 . l-s 2 y~/ i 2 \ / i 2 V „| 



Now being based on the result (|A4j) we introduce the function 



F w (2,y) 



5 



1 + + 7:2/ 



1 2 9 2 



<9y 6 dy 2 



2v Jl-2v 



-1 



^2 _ 2 ^ + ^ 



whose discontinuity at z = x reads 

= d isCi « (i + , | + i, 2 ^) /_; ^ | - , 2 )p, g) 
-H)(^)EWl'(;) ffi (;). 



The asterisk in the sum in ( |A6| ) denotes that for z/ = the terms with Z = and I 
actually absent. 

6 The discontinuity of the function f(z) is defined as disc 2=2; f(z) = (f(x — iO) — f(x + iO)) = 
iO) 
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Eq. (|A6|) is the natural building block for the desired convolution kernel. Indeed, accord- 
ing to the definition of the gluon forward like functions G 2u (y, t) with n = 1v + I — 1 



^2i_ 2 „(*)= f 1 dVV n G% ) (!J.t) witl. n>2u-l. odd 
Jo 



(E)< 

y,vj vviuii io^*v — ±, u<au.. 

(A7) 



Using ( |A7| ) together with ( |A6| ) it is straightforward to check that the integral convolution 

oo „i 

V / dy-ldiac^F^^y) +disc, = _,F^(z, 2/ )}G£ ) ( 2 /,t) (A8) 

results in the formal series fllSf) - 

The trick that allows to derive the expression for the convolution kernel consists in the 
explicit calculation of the discontinuity of F^ 2u \z,y) at z = x stemming from the cut at 
1 ~~ \A ~ £ 2 < x s < 1 + a/1 — ^ 2 and from possible poles at z s = 0: 



f dyG%\y,t) disc z=x F^(z,y) 
Jo 

d_ 

dy"' 6" dy 2 J J_^" ^-xl-i 



(A9) 



Finally, the expression for the gluon GPD H 9 ^ through the set of the forward like 
functions reads 

oo r \ 



H^ E \x,^t) = / dy-{disc z=x F^\z,y) +d\sc z= . x F^\z, y)} G { 2 E J(y,t) 

u=0 ^° 

= E V [ R9 {E) M i*>t>t) + H" ^ M(-x,Z,t)] 



1-72 ) &l_ 2 [-)B^> (t) 



(A10) 

where the functions H 9 ^ E ^ v \x,^,t) defined for — £ < x < 1 are given by the following 
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integral transformations: 



(s>0- / d v 



71 



1/ 9 1 J 2 

3 l 1 -^ + 2 2/ ^ 



<£)- / dy 
^ Jo 



1 / 9 1 J 2 
3 [ 1 '% + 2 y W 2 



, 9 x2 2^-3 



6y 



82 x 2 - 2u (l-s 2 ) 



f 2x s -x 2 s -e 

83 x 2 ' 2u (l - s 2 ) 



si 

2i^-i-2 



'2x s -xj-g 2 



2,-1-3 \ t)t p i( (2v-l)(2v-l + l) 



(All) 



with P- n (x) = P n -i(x)- Note, that in ( |A11| ) we employ the standard notations adopted 
in [T5|, P^] . Namely, x s = 2 (jz^fe , Sj, (z = 1, ... 4) stand for the four roots of the equation 
2x s — x 2 — £ 2 = given by the following expressions: 



Sl 
S3 



y.-y/{\-xy) (l + /i 2 )-(l-?/ 2 )); 
fj,+ y/(l-xy) (l + /i 2 )-(l-y 2 )); 



- (A - y/(l-xy) (l + \ 2 )-(l-y 2 )) ; 



- A + V(l-xy) (1 + A 2 ) - (1 - y 2 ) 



(A12) 



where 



(A13) 



yo and ^ are the two solutions of the equation s x = s 2 ; while yi and i are the two solutions 
of the equation S3 = S4; 



Vo 



x(l + /i 2 ) /x 2 (l+/i 2 ) 



2^2 



H 2 ; 



(A14) 



yi 



x (1 + A 2 



, x 2 (l + A 2 ) i 



-A 2 . 



(A15) 
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B. POLARIZED GLUON GPDS 



As usual, to sum up the formal series for H 9 (^0]) we introduce the set of polarized 
gluon forward like functions AG2 V (y,t) whose Mellin moments generate the generalized 
form factors B^^t). with n = 2v + I — 1 

-i 



dyy n AG2u(y,t) with n > 2, even. 



(B16) 



The resulting expression for H 9 through AG2u(y,t) reads 



u=0 



<9z dU 2 



(B17) 



where H 9 ^(x,£,t) defined for — £ < a; < 1 is given by 



= 9(x > 0- / dy 



9 1 2 d 2 
-A l ~ V dy- + 2 V W 



AG 2v {y,t) 



8 2 



ds 



x 



2-2u, 



1-s 2 



+6(\x\ < £)- / dy 



71 







K i_!/ l + ^ |AG2 - to ' t) 



a/2x s - - £ 2 



ds 



i V 2 ^ - x l - £ 2 



£ cl '-"(l) w (? 



Qy 



2v-l-2 



(2u-l)(2v-l + \) 



(B18) 



For even N the polynomiality condition (|ToD require 

r 1 N 
/ dxx N - 1 H 9 (x,^t) = Y,t k h 9 Nk {t) = 

JO fc=0 



2V n+1 



f N EE^ ^ » (! + ") (2 + n) (3 + n) 1(f) I \ 



n = 2 i=l 
even odd 



C; 9 ■ 2^ r(l + =*±Z) r(| + =gifM ) 



(B19) 



The corresponding set of coefficients h 9 N k (t) is expressed through the generalized form factors 
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B 9 nl {t) as follows 



JV n+1 



fc+i-iV+l 
2 



•i + fe — jv) r 



2-fc+l+jV' 



n = 2 1=1 
even odd 



3 . 2 fc+ir( 1 + fc +'- jV ) r(2 — k + N) 



X 



n (1 + n) (2 + n) (3 + n) r(7V) 

r(^v) r(5±S) 



(B20) 



The expression for the polarized forward like function AGo through the t-dependent 
polarized gluon density Ag(y,t) reads 

AG (x, t) = - 9 -x 2 [ %Ag{y, t) + 3x [ %Ag(y, t) + \ / ^A^(y, t) . (B21) 

^ J x y J x y ^ J x y 



Finally, to sum up the formal series for H 9 ( ps ^ (^l|) the set of polarized gluon forward like 
functions AG^f\y^t) whose Mellin moments generate the generalized form factors B^^t). 
with n = 2v + I 

• i 



B 



9 (PS) 



nn-2u(t)= I dyy n AG™(y } t) with n > 2, even, 
o 



(B22) 



The resulting expression for H 9 through AG^ reads 

00 £2k 



H 9 (x, £, t) = V &9 iPS) ^)-H 9 {PS) H t) 



u=0 



where H 9 ^ PS ^ u \x,^,t) defined for — £ < x < 1 is given by 



(x > £)- / % 



7T 



+0(|s| <£)- / 



7T 







AG\ U J {y,t) 



*2 



(is 



1-2^/ 



(B23) 



1-^ 



3 V ^ dy 2^ dy 2 



s-.i 



^J2x s -x 2 -£ 2 
xl- 2v {l-s 2 ) 



ds 



, y/2x s -x 2 -e 



(n N 2 2I/-2 



-/-2 it!'' 



6y 



£/ (2i/-Z + 1)(2i/-Z + 2) 



(B24) 
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For even N the polynomiality condition (|10|) require that 



N 



dxx N ^H^ s \x^t) = Y,eh a N {p k s \t) 

even 

,,fv &>™mp f 1 } n(1 + n) {2 + n) (3 + n) r(i)r(A ° 



n=2 Z=0 \ > / ^ ^ ""-V- 1 - 1 2 / V 2 

even even 

The coefficients h 9 j!f k if) at powers of £ in (|B25| ) are given by 

19{PS) (A _V^Y^ 5s(PBf)/^/_ 1 \4±ipK M 2 j 

'^iv.fe w 2.^ 2-^ D m L > o 2^ 2 r( 2+fc +'~ jV ) r(i fc i A r ) 

even even 

n (1 + n) (2 + n) (3 + n) r(iV) 



(B26) 



The expression for the pseudoscalar forward like function AG[ ) P ' S ' ) through the pseudoscalar 
combination of t-dependent polarized gluon densities g( ps \y, t) = Ag(y, t)+rAe 9 (y, t) reads 

AGi PS \x,t) = ^ f%Ag^Xy,t)-S X f %A g ^\y,t) -§ f" ^Ag^ s \y,t) . 

^ j x y J x y ^ J x y 

(B27) 
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